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Abstract. We study the spherical collapse in the Parametrized Post-Friedmannian (PPF) 
scheme. We use a general form of the PPF parameter related to the Poisson equation and 
found the equations to solve that includes a non-trivial fifth force coming from the convolution 
of the modified gravity term in the fc-space. In order to compute a concrete model, we use 
the parametrization proposed by Bertschinger and Zukin. The equations of the spherical 
collapse are solved assuming a Gaussian density profile and we show there is no shell crossing 
before reaching the turn around point. We show that the fifth force does not satisfy the 
Birkhoff’s theorem and introduces different behaviors for the density threshold 6 C , which in 
this case depends on the size and shape of the initial density profile, and therefore one expects 
a different statistic of the collapsed objects in the universe. 
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1 Introduction 

One of the most exciting open problems in cosmology is the well understanding of the mech¬ 
anism that carries the observed late-time accelerating expansion of the Universe [1, 2], Since 
the first evidence of the need of a new component in the universe called Dark Energy (DE), 
the introduction of the well known cosmological constant (A) as DE has succeed a variety of 
observations [3] and has established the standard cosmological model ACDM, which assumes 
General Relativity (GR) and also introduces a dark matter component [4]. However, the 
cosmological constant leaves doubts about its nature. Hence, many other theoretical alter¬ 
natives have been also proposed that try to explain the accelerating expansion in different 
ways. One approach is to add an “exotic matter” source with negative presure in the right 
hand side of the Einstein equation; within this approach one has quintessence [5], k-essence 
[6], and perfect fluid models [7]. The other approach is the called “modified gravity” (MG) 
models, that modify the left hand side of the Einstein equations. This includes / ( R ) gravity 
[8-10], scalar-tensor theories [11-13] and braneworld models [14, 15]. 

If any of the proposed theories attempt to be valid and to shed light for a better under¬ 
standing of the dark sector then they must be confronted with observations. Typically, one 
has to work out the perturbation equations for each theory and incorporate them into an 
Einstein-Boltzmann solver or a N-body code to construct a list of observables to compare to 
the data. This task can be cumbersome and time-consuming. Over the last years it has been 
of great interest the developing of an intermediate step between theory and data that allows 
us to avoid working out each model and encompass a wide range of theoretical schemes, at 
least at background and first order perturbation level. This phenomenological approach has 
been dubbed as “Parametrized Post-Friedmannian” scheme (PPF) [16], inspired by the so 
called Parametrized Post-Newtonian (PPN) formalism [17] created to confront alternative 
theories of gravity with Solar System measurements. The idea is to assume a ACDM model 
at background level and to parametrize deviations of GR at first order in perturbations. 

There are several modified gravity parametrizations, of which one widely used is the one 
proposed by Bertschinger and Zukin [18] (BZ), where two functions of time and scale are 
introduced motivated by / ( R ) solutions in the quasi-static limit which has been shown in 
recent works to be valid for large scale structure applications in / ( R ) models [19] . This is 
the approach we closely used in this work. In addition to this parametrization some other 
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general frameworks have been developed, such as the Hu-Sawicki model [16], its relativis¬ 
tic corrections [20] and the quite general approach proposed by Baker and Ferreira in [21]. 
These parametrizations are not independent, but related to each other as mapped en [22], 
The target that all these different approaches have in common is to create a link between 
modified gravity theories with its corresponding PPF parameters. 

However, the objects from which cosmological information is extracted can be beyond the 
linear perturbation regime (galaxies and clusters of galaxies). Therefore, the study of grav¬ 
itational collapse and dynamics of collisionless systems at non-linear perturbations is neces¬ 
sary to complete the description of structure formation in the Universe. Although, the full 
non-linear dynamics has to be dealt only with numerical simulations, one can make some 
assumptions about the system symmetry in order to construct some analytical models that 
provide a valuable first approximation to the non-linear regime. One of these analytical 
models is the so called Spherical Collapse (SC) model introduced in [23]. In this model a 
spherical, symmetric overdensity evolves to form a virialized bound system via gravitational 
instability. This model can be applied in different cosmological contexts to study the abun¬ 
dances of virialized objects through the Press-Schechter formalism [24, 25] or relaxing the 
spherical symmetry through its extension of Sheth and Tormen [26]. 

Lately, many authors have used this model to study the halo mass function and the linear 
halo bias [26] in different cosmological scenarios. Adding the presence of massive neutrinos 
[27], in braneworld cosmology [28], in Galileon gravity [29], in context of f (R) theories [30- 
32] or more recently in chameleon theories [33]. In this work we construct and solve the SC 
model in the BZ scheme developed in [18], but after constructing the formalism one can move 
on and solve the SC in any other kind of parametrization, which is the main reason of using 
a parametrized scheme of modified gravity. 

This paper is organized as follows. In Section II we present the PPF scheme and the / ( R ) 
case is discussed in Section III. In Section IV we implement a PPF scheme to the SC Model 
and we show our main results. Finally, we summarize our findings in Section V. 

2 Parametrized Modified Gravity 

Considering linear scalar perturbation and the perturbed metric tensor in the conformal 
Newtonian gauge: 


ds 2 = a (ry) 2 [- (1 + 2T) dry 2 + (1 + 24>) dx 2 ] , (2.1) 

where ry refers to the conformal time. The growth of structure in the universe at linear 
perturbation level is determined by the Poisson and the gravitational slip equations. These 
two equations in Fourier space read (neglecting any contribution from the anisotropy stress, 
a): 

k 2 m = —AirGa 2 pA (2.2) 

<£ + T = 0. (2.3) 

where k is the wavenumber, A = 5 + 3aH(l + uj)6/k 2 is the comoving density contrast and 
a is the scale factor. Note that equation (2.2) is not the original Poisson equation for the 
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potential $ that comes from combining the time-time and the time-space components of 
Einstein equations. Instead, we have derived an analogous Poisson equation for the potential 
'I', which is more closely related to observation. In fact, T is the potential that governs the 
non-relativistic particles motion. 

Alternative theories of gravity generate a different cluster dynamics and CMB evolution 
from the ACDM model. Equations (2.2) and (2.3) have in general different forms, depending 
on the class of Modified Gravity (MG) theory involved. Thus, the task to study the evolution 
of structure and CMB anisotropies in different MG models could be tedious and cumbersome. 
For this reason the PPF formalism was put forward that let us study different theoretical 
MG models within the same first order perturbation formalism. This can be achieved by 
introducing two free functions of time (or scale factor) and scale: p(a,k) and 7 (a, fc), that 
modify the Poisson and gravitational slip equations, in order to encompass a group of MG 
models. Accordingly, one can write: 

/c 2 'P = —AnGcPp (a, k) pA (2.4) 

<P + 7 (a, k) T = 0. (2-5) 

where the usual GR case is recovered when p = 7 = 1. A general modification of gravity 
through an arbitrariness in the explicit form of the functions p (a, k) and 7 (a, k) could be 
tedious and difficult to constrain. Silvestri et. al. [34] have shown that in a set of viable 
models in which local theories of gravity are considered possessing only one extra degree of 
freedom and second order equations of motion in the Quasi-Static Limit (QSL), the functions 
p (a, k) and 7 (a, k) must be ratios of even polynomials of second order in k and in general the 
numerator of p(a,k ) is the same as the denominator of 7 (a,k). Introducing five functions 
of the background pi (a) one has: 


7 = 


P = 


Pi (a) + P 2 ( a ) k 2 

1+P3 (a) k 2 

( 2 . 6 ) 

l+p 3 (a) k 2 

Pa (a) + p 5 (a) k 2 ' 

(2.7) 


Although this form of parametrized MG is found in the QSL, in which k/aH 1 is assumed, 
these expressions are quite suitable for near-horizon scales since near and super-horizon effects 
have practically no impact on observables in viable MG models when relativistic effects are 
neglected [35] ; for a treatment of relativistic effects see [20] . 


3 Parametrized MG in / (R) Theories 

One of the simplest MG models is the / ( R ) gravity in which the 4-dimensional action is 
given by some general function / ( R ) of the Ricci scalar R, in the Jordan frame: 

s = 7^2 / ^ x V~9f ( R ) + S rn {Pfiu, Ipm) , ( 3 - 1 ) 

where k 2 = 8irG, and S m is a matter action with minimal coupled matter fields ip m . By 
varying the action (3.1) with respect the metric one obtains the field equation: 
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IrR^u - (R) - V pSvfR + g^Ofn = k 2 T^, (3.2) 

where /r = df/dR and T /JL/ is the matter energy-momentum tensor. The trace of the above 
equation is given by 

3D/i? + /rR — 2/ (R) = k 2 T, (3.3) 

where T = g^ v T^ v . Equation (3.3) governs the dynamics of a scalar field <p = fa, dubbed 
scalaron [36], which is a propagating scalar degree of freedom. Any viable / ( R ) model must 
satisfy a number of conditions, i.e., to avoid instabilities and tachyonic scalaron (Jrr > 0); 
to evade the existence of ghosts [37] (/r > 0); to recover GR at early times (/ (R) —> R — 2A 
for R —> oo) and at solar system (|/j?o| — 1 < 10~ 6 , /r evaluated today [16]. See [38] for a 
complete discussion about / ( R ) theories. 

Given an expansion history H (a) and H (a), it is possible to design a / ( R ) that solves 
(3.2) constrained by the boundary condition /rq. Thus, / ( R ) models constructed in this way 
can have an expansion history practically indistinguishable from the ACDM model [16, 39]. 
To be able to distinguish among MG models, it is then necessary to study the observables 
of linear and non-linear theory, e.g. the evolution of cosmic structures or the growth factor, 
preferably at different redshifts, and to compare them with observations. 

To describe structure formation in a / (R) model it is necessary to expand Jr in per¬ 
turbation theory, 6/r. Hence, in / ( R ) new terms appear in the standard Poisson and 
gravitational slip equations. One obtains, at first order [35], 


k 2 ^! 


>SfR 

2 Ir 




h 2 ) s -R 

' JR 




Ir 

Ir 


T + 


a 2 p 
2 k 2 Ir 


A, 


SfR 

Ir 


(3.4) 

(3.5) 


As discussed in [40], in a 
27 t /mf R , where 


/ ( R ) model the Compton wavelength of the scalaron Ac = 




(3.6) 


is the effective scalaron mass, that sets two regimes in the dynamics. In scales above A c 
modifications are negligible and GR is recovered. On the other hand, in scales smaller than 
the Compton wavelength the growth is enhanced by a fifth force and in addition the metric 
potentials are no longer equal. Also, if one assumes the QSL, the two above equations can 
be reproduced by equations (2.4) and (2.5) with a suitable p and 7 : 




1 1 +(4/3 )Q 

Ir 1 + Q 
1 +(2/3 )Q 
1 + (4/3) Q ’ 


(3.7) 

(3.8) 


in which Q = Inspire by this in [18] it is proposed the following parametrization for 

p and 7 : 
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(3.9) 




BZ 


7 


BZ 


1 + §A jk 2 a s 
1 + A \k 2 a s ’ 
1 + | A \k 2 a s 
1 + |A \k 2 a s ’ 


(3.10) 


where Ai represents the Compton wavelength and s encodes its time dependence. In [35] it is 
shown that this last form of parametrizing with s = 4 is quite accurate and can be safety used 
for deriving constraints on / ( R ) models for next generation of large scale surveys. Notice 
that (3.9) and (3.10) fall within the more general form of p and 7 discussed before, eqs. (2.6) 
and (2.7). 


4 Spherical Collapse within the PPF formalism 

In this section we will develop a semi-analytic SC model considering a MG force due the 
modified Poisson equation (2.4) where p (a, k ) will be given by equation (2.7). The SC model 
treats the evolution of thin, adjacent concentric mass shells forming a constant spherical sym¬ 
metric perturbation in the cosmic density field. See [41] for a standard treatment of SC model. 

In the Newtonian limit the evolution of each shell is governed by the Newtonian force law. 
Considering dark matter and dark energy components one has: 

= _VT - ^H 2 n m r + H 2 n A r, (4.1) 

where r is the physical position of inner shells, f l m and are the dimensionless matter 
and cosmological constant energy densities, respectively. Usually, the SC model considers a 
constant, or top-hat, initial density profile, which is saved from shell crossing in both Einstein 
de-Sitter (EdS) and ACDM model. Similar as it happens in / ( R ) models, an initial top-hat 
profile suffers from shell-crossing generating a large spike near the perturbation edge [31]. To 
avoid shell crossing, in general, a perturbation with a monotonically decreasing Si profile is 
needed: 


S (r) = f S (/) r' 2 dr'. (4.2) 

r 6 Jo 

where the bar denotes volumen average. For this reason we will consider an initial Gaussian 
profile for the density profile. This alleviates the shell crossing and incorporates a more 
realistic perturbation, 


Pi (r) = po exp (—A 2 r 2 ) , (4.3) 

where po is the initial magnitude of the density profile at r = 0. Before proceeding to solve 
(4.1) one has to calculate the explicit form of the force term —V'k, which is the only place 
a MG model is captured, assuming a standard ACDM background dynamics as in the PPF 
formalism; otherwise H also encodes a different background dynamics. 

Firstly, we have to transform the Poisson equation (2.4) into the configuration space, in the 
QSL we have: 
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V 2 T = AnGci 2 (dp * F~ l p) , (4.4) 

J 7-1 means the inverse Fourier transform and denotes the convolution between the density 
contrast and the p function. Assuming (2.7), we obtain: 


F p = p(r) = 


1 Pi ( 1 Vi\ e 


4vrp 5 \p4 p 5 


1 


= -/3r 


= -T-POL — 
47 t r 


+ A5 d (r), 


where (5 b denotes de Dirac delta and 


a = — - — 

~ Pi Pb 

P5 

A - V P 3 (°) 

A — Inn p = ——. 

k— KX) P5 (a) 


+ A5 d (r) 


Expanding the convolution in the Poisson equation, we finally get 


(4.5) 


(4.6) 

(4.7) 

(4.8) 


V 2 tf = 4 t tG 




dr' + Ap — (a + A) p 


(4.9) 


where we have defined x( r i r ') = (| r — r/ |) sin (0) F (r) = /3 2 a^r— and M (r) 

is the mass function that provides the total mass inside a radius r. Integrating this last 
expression to obtain the wanted force expression: 


- V’P = -- 


4vrG 






+ Hfalr. (4.10) 


Setting a = 0 and A = 1 causes the first term of the right hand side (r.h.s) to vanish and 
the usual force is recovered. This modified force not only implements a time dependence 
transition between standard gravity at larger scales and A-times stronger gravity at smaller 
scales, but also introduces an important phenomenological effect that comes from the con¬ 
volution operator in the Poisson equation. This operator triggers off the mixing of scales 
at some physical radius i.e., the force at some physical position r not only depends on the 
matter inside that radius but also on the configuration of the matter outside. Although, this 
scale dependence is weak because of the p function is quite cusped at the origin, this breaks 
the Birkhoff’s theorem and introduces a series of effects that alter the way an spherical pulse 
collapses, and therefore, a different value for the density threshold 6 C is now obtained. The 
main reason behind this is that the collapse now depends on the configuration (size and 
shape) of the density profile due to the scale dependent new force. Given this, an initial 
homogeneous matter density evolves to a non-constant profile that eventually may suffer 
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Figure 1. Evolution of the innermost shell in different gravity models. In blue (solid) the ACDM 
model with standard gravity; in magenta (dashed) and in green (dotted) MG gravity models with 
a = —0.001 and a = —0.0025, respectively, keeping /3 q = 1 Gpc -1 in both models. 


from shell-crossing. This the reason why we assumed a decreasing density profile (4.3), for 
which we have to follow the evolution of each concentric shell in order to avoid shell-crossing 
before the turn-around time. A similar situation is studied in [42], where a Yukawa force is 
considered. 


Now we proceed to solve numerically (4.1) to find the radius of each inner shell as 
function of the scale factor, r(a). To accomplish this, it is convenient to make a change 
variable r = —x, as in [43], in which r, is the initial position of the corresponding shell and 
a* is the value on the scale factor at initial time, that we set to be a* = 0.01 (z ~ 100) 
that is well inside matter-dominant era and plenty before the accelerated cosmic expansion 
(z < 1). Also, we set the perturbation amplitude ~ 10 —2 , to have initial conditions within 
the linear theory. Taking into account the results of the previous section, the equation of 
motion becomes 


cl 2 x 
d t 2 


47 tG a? 


■i ri 


rXi 


X {x'\ x') 



x" 2 dx" 


Integrating by parts the first term of the r.h.s one gets: 


1 H 2 Q° 

1 a n 0 iL m 

2 x 2 


(l + + HqXI^x. 


(4.11) 


d 2 x 

dt 2 


2irG a 
x 2 rf L 


3 px 

i M (Xj) / xix.X^x'^dx'- 

i L 40 


i et2q0 

(i + 3j) + h 2 q° a x, 



x" 2 dx"dx 


(4.12) 


where we defined Xi = ^ B-i ■ 

' i 
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Figure 2. Scale factor at turn around as function of the initial shell position, showing that there is 
no shell crossing, at least not before the turn around time. In blue (solid) ACDM model; in magenta 
(dashed) and in green (dotted) MG gravity models with a = —0.001 and a = —0.0025, respectively, 
keeping /3 q = 1 Gpc -1 in both models. 


To perform the integration we have to specify the MG model, given by (2.7). We used 
the BZ model, eq. (3.9), to determine (4.5), in which 

a = l-A, (4.13) 

P 2 = Pla~ A = A rV 4 , (4.14) 

in which we have taken the values of A that imply small deviations from GR and we set 
Ai = 1 Gpc. We have assumed a background ACDM model with no curvature and = 0.31. 
In the figure (1) we show x (a) for the innermost shell for ACDM model (a = 0) and for two 
particular MG models with a = —0.001 and a = —0.0025. As it can be seen from the figure 
each shell with initial radius r* expands until it reaches a maximum and then turns around 
to eventually collapse as r —>• 0. For a < 0 the gravity is stronger than the ACDM model 
causing the shells to collapse faster and for a > 0 a slower collapse takes place. One refers 
as the turn around time ( t ta) to the time at which the outermost shell reach its maximum 
expansion. After turn around the collapse follows but, physically, the perturbation will never 
reach the singularity due to merging or violent relaxation that establish a virial equilibrium. 
Therefore, it is customary to assume that the density contrast of the dark matter halo virial- 
izes at twice the turn-around time. From figure (1) it is clear the known effect of a fifth force 
to increase the gravitational potential to develop an earlier collapse than the ACDM model. 

Because of a non-constant initial density profile the evolution equation for the inner 
shells are no longer self-similar. For instance, in the top-hat case each shell evolution equa¬ 
tion is indeed self-similar, meaning that each inner shell reaches its turn-around point (and 
its collapse point) at the same time. In the case of a monotonically decreasing initial density 
profile, as figure (2) shows, the time at which each inner shell reaches its turn-around point 
(and its collapse point) is a monotonically increasing function of the initial shell position r % . 
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For the sake of not spoiling the SC model we need to set the edge of the initial Gaussian 
profile so that the outermost shell reaches the turn-around point before the innermost col¬ 
lapses. In general, this can be achieved demanding Si (0) /Si ( Ri ) < 1.587 [43], where Ri is 
the radius of the initial density contrast. In the case of a Gaussian profile this condition can 
be easily accomplished by setting the edge of the Gaussian pulse at Ri = 0.8A" 1 , see eq. (4.3). 

Figure (3) shows the evolution of density contrasts a few time steps before the inner¬ 
most shell collapses; it is the latest time the SC model is still valid. The top-left panel 
shows the ACDM model, the top-right and bottom-left panels show MG gravity models with 
a = —0.001 and a = —0.0025, respectively, keeping /3q = 1 GpW 1 in both models, see eq. 
4.14. As it can be notice from the figure the three models (ct = 0,-0.001,-0.0025) have a 
similar behavior with the only difference in the collapsed times. In the MG cases the collapse 
is slightly faster than in ACDM model (or slower if a > 0). This can be seen in the bottom- 
right panel where the three models are plotted at the same time (the time at which the 
innermost shell of the model with a = —0.0025 collapses). As expected, the gravity model 
with a stronger deviation from GR collapses faster, thus, 5 is bigger. Despite the differences 
among models, when they are plotted at their collapsed time (the time at which their in¬ 
nermost shell collapses) the final behaviors look more similar, finding the same pattern: a 
model with stronger deviation has a denser profile, according to the a parameter. But the 
small differences shown in figure (4) cause considerable changes in the density threshold S c , 
as it is shown below. 

The final profiles we have obtained do not represent virialized structures, because their 
outermost shells have just past the turn around time and they have not reached a virial 
condition. Following, we test how much the collapsed structures differ from a typical collapse 
profile. Accordingly, we fit a Navarro-Frenk-White (NFW) profile [44] to our final profile, 
which is given by: 


P = 


Po 



(4.15) 


where c is called the concentration parameter and r v is the radius of the virialized halo, and 
their values vary from halo to halo; po is a typical halo density that depends on r v , c, and 
cosmological parameters. Figure (5) shows the numerical density profiles obtained (dots), 
divided by the critical density evaluated at a = 1; and the NFW fittings (solid). As it can be 
seen from the figure, right panel, the fittings are very close to the density profiles (obtained 
from the SC model) in the innermost shells, but they start to deviate from r ~ IMpc. This is 
because the outermost shells of each pulse are in their way to collapse and have not reached 
a virial condition yet. 

Furthermore the collapse time, or the time at which the halo virializes, depends on the 
initial value of the density contrast Sq] the higher 5q the earlier the overdense pulse collapses. 
The value of the linear density contrast extrapolated at the time of collapse is referred as 
the critical threshold S c . This quantity is important because it represents a key element in 
calculating the mass function [24], the number density dn of collapse halos with mass in the 
d M range. The analytic expression is given for S c in the case of the EdS model [41]: 
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Figure 3. The latest time steps of the density contrast 5 evolution before it reaches the collapse 
criteria for the innermost shell. Top-left: ACDM model; top-right: MG model with a = —0.001 and 
/?o = 1 Gpc” 1 ; bottom-left: MG model with a = —0.0025 and /3q = 1 Gpc” 1 ; bottom-right: The 
previous three models at the same time (the time at which the model with a = —0.0025 collapses). 
In blue (solid) the ACDM model, in magenta (dashed) and in green (dotted) the MG models with 
a = —0.001 and a = —0.0025, respectively. 


4 Oc) = ^vr 2/3 [H (a c ) t c ] 2/3 , 

and for the ACDM model 


(4.16) 


S c iflc) 






1 + Si 



(4.17) 


g (a c ) is the linear growth factor evaluated at collapse time. In MG models one can use the 
same expression for 5 C as equation (4.17), but in this case g (a) is the linear growth factor for 
the modified linear differential equation for the evolution of 5 in the PPF formalism, which 
in Fourier space is given by 


5" + US' - (a, k) 5 = 0, (4.18) 

where ' denotes derivation with respect conformal time. This last equation follows from 
the energy-momentum tensor conservation, considering MG only enter through the Poisson’s 
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Figure 4. The density contrasts <5 at their collapse time, given the corresponding criteria for the 
innermost shell. In blue (solid) the ACDM model, in magenta (dashed) and in green (dotted) MG 
models with a = —0.001 and a = —0.0025, respectively. 



Figure 5. The density profiles (over the critical density today). Left panel: in blue the ACDM 
model, in magenta and in green MG models with a = —0.001 and a = —0.0025, respectively. Dots 
represent density profiles obtained from the SC model at their collapse criteria for the innermost shell 
and in solid lines NFW fitting profiles with r v = 3.19 and c = 4.95; r v = 2.99 and c = 6.52; r v = 2.86 
and c = 8.60, respectively. Right panel: the same as in left panel but in log scale. 


equation. 

Using the BZ parametrized MG model (3.9), Ai is the only free parameter which can be 
related to the / ( R ) dimensionless parameter B$ [45]: 


B 0 = 


Ir 

fRR 


ml \ 2 _ 0-2 

^r Al - -#-Po 


(4.19) 
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Figure 6. Linear density threshold at collapsing time, S c , for different models. In dotted (blue) the 
EdS model, in solid (blue) the ACDM model and in dashed lines (green) MG models with different 
a values, from bottom to the top: -0.001, -0.0025 and -0.005. As a decreases its values, 5 C departs 
more from the ACDM model. We set /3q = 1 Gpc -1 in all MG models. 


In [46] the authors obtained an upper limit to this dimensionless length scale Bq < 0.14 at 
95% c.l. from CMB measurements. From this bound we have that (3o > 1 Gpc -1 , which 
means that in scales larger than Ai ~ 1 Gpc one has to recover GR. 

In figure (6) results are shown for 5 C ( a c ) varying the parameter a. The usual constant 
result (1.686) for an EdS Universe is shown in dotted line (blue) and in solid line (blue) for 
the ACDM model in accordance with known results [47]. For the MG models presented here 
we note that S c (a c ) tends to the ACDM case as a —>• 0. As a o decreases (starting from 
zero) S c (a c ) deviates from the ACDM case and for much larger values it gets over its EdS 
counterpart. 


5 Conclusions 

We have worked out the SC in MG models of the PPF type. We found some analytic ex¬ 
pressions for a general modification of gravity of type given by equation (2.4) and found 
numerical solutions for the parametrization designed to mimic the / (R) gravity behaviour 
proposed by Bertschinger and Zukin (BZ). As it happens in / ( R ) models, an initial constant 
profile suffers from shell-crossing, making our analysis invalid. Therefore, in order to avoid 
shell-crossing a Gaussian density profile was proposed, but in general a decreasing profile 
works as well; figure (2) shows shell-crossing is avoided in the models presented in this work. 
We solved with some detail the evolution of two different MG models with a = —0.001 and 
a = —0.0025, respectively, keeping the MG transition scale fixed to (3q = A] -1 = 1 Gpc -1 . In 
figure (3) and in figure (4) show, as it was expected, gravity models with strongest deviation 
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from GR collapses most and fastest. 


Though the final density profiles obtained do not represent virialized structures, we 
fitted a NFW profile to the each of them. In figure (5) we showed that the NFW fittings are 
very close to the final density profiles (obtained from the SC model) in their innermost shells, 
but they start to deviate for the most external shells. This is because the outermost part of 
the pulses have just passed the turn around point and they are in the way to reached a virial 
condition. It is interesting that NFW profiles perform good fits to MG collapse models, at 
least at the stage considered in this work. 

Finally, figure (6) shows the density threshold 5 C for different MG models. These can de¬ 
viate considerably from the ACDM behaviour depending upon the values of the a parameter, 
which dictates the strength of the fifth force. We kept the value of /3 q fixed to /3q = 1 GpcU 1 
which determines the time transition between standard gravity force at early times and A- 
times stronger gravity at recent times. Further, we noted that the fifth force introduced does 
not follow the Birkhoff’s theorem making the spherical collapse model to depend also on the 
density profile configuration. Our results indicate that the collapsed objects are denser than 
ACDM predictions and this effect strengthens at smaller redshifts. 
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